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Abstract 

An efficient technique to solve precision problems consists in using exact computations. For 
geometric predicates, using systematically expensive exact computations can be avoided by 
the use of filters. The predicate is first evaluated using rounding computations, and an error 
estimation gives a certificate of the validity of the result. In this note, we studies the statistical 
efficiency of filters for cosphericity predicate with an assumption of regular distribution of the 
points. We prove that the expected value of the polynomial corresponding to the in sphere 
test is greater than e with probability 0(elog -) improving the results of a previous paper 
Keywords: Computational geometry, Delaunay triangulation, exact arithmetic. 



1 Introduction 

The assumption of real-number arithmetic, which is at the basis of conventional geometric 
algorithms, has been seriously challenged in recent years, since digital computers do not ex- 
hibit such capability. Geometric algorithms involve the evaluation of predicates; to guarantee 
the structural correctness of the results, predicates must be evaluated exactly. A geometric 
predicate usually consists of evaluating the sign of some algebraic expression. In most cases, 
rounded computations yield a reliable result, but sometimes rounded arithmetic introduces 
errors which may invalidate the algorithms. Assuming error-free input data, the rounded 
arithmetic may produce an incorrect result only if the exact absolute value of the algebraic 
expression is smaller than some (small) e, which represents the largest error that may arise in 
the evaluation of the expression. The threshold e depends on the structure of the expression 
and on the adopted computer arithmetic. This is basically the philosophy behind the notion of 
arithmetic filters, whose function is to adjust the arithmetic overhead, so that no more effort 
is expended than required by the test instance. 

It is therefore of interest to estimate the frequency with which recourse to arithmetic 
engines more powerful than standard platforms is necessary. Such analysis must be carried 
out by making some a priori hypothesis on the distribution of the input data, which are 
treated like random variables. Since for our objectives only the absolute value of the algebraic 
expressions is significant, hereafter "value" is to be intended as "absolute value". 
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In a previous paper [DP98 , we have carried out such analysis for two crucial geometric pred- 
icates, the orientation test (which-side of a hyperplane) and the insphere test (inside/ouside a 
hypersphere), on the hypotheses that the input points were uniformly distributed either in the 
unit ball Bs or in the unit cube Cs = [— 1, 1]* in <5- dimensional space. Our results were that, 
for a small value V, the probability that the result of the orientation test is < V is 0(V) in 
all dimensions, whereas for the more complex insphere test we obtained bounds sublinear in 
V . Specifically, we obtained 0(V 2 ^ 3 ) in dimension 1 (which is tight), 0(V rl/ ' 2 ) in dimension 
2, and 0( V 1 / 2 In V) in higher dimension. 

Later on , we discovered a discrepancy between these theoretical findings for 5 > 1 and 
the results of extensive simulations, which seemed to exhibit a linear behavior (see below). 
This observation motivated a finer analysis, reported in this note, whose conclusion is that for 
5 > 1 and for 8 + 2 points pi,ps, ■ ■ ■ ,Ps+2 uniformly chosen in the unit ball, the probability 
that the value of the determinant, embodying the insphere test of ps+2 versus pi,pa, ■ ■ ■ ,ps+i, 
is < V is 0(Vln(l/V), in closer agreement with the simulations. The results extend to 
points uniformly chosen in a cube. We also present an application of this analysis to the 
three-dimensional insphere test carried out with floating point arithmetic. 



2 Analysis of the insphere test 

The algebraic expression embodying the predicate which tests if a point ps+i belongs to th e 
sphere S passing through points P1P2 • • -Ps and the origin, is the following determinant [DP96|: 
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As mentioned in the Introduction, in dimension 1 the insphere test reduces to an in-interval 
test and is only of moderate interest. Nevertheless, we have obtained the following tight bound 
jDP98t Pro6(|Ai| < V) < ^^Vi ~ 5.36V 2/3 

We now turn our attention to higher dimension, and let c = (^-, . . . , denote the center 
of the sphere S. In the above determinant, subtracting column i times a from the last column, 
enables us to rewrite A^ as 
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This expression can be more synthetically rewritten as W — \cps+i \ 2 —r 2 , i.e., W is power(ps+i, S) 
of point ps+i with respect to the sphere S. Notice that power (ps+ 1, S) is positive if ps+i is 
external to S and negative if it's internal. Therefore random variable As is the product of the 
two random variables \p1p2 ■ ■ - Ps\ and power(ps+i, S) ( of which, incidentally, \p1p2 ■ ■ ■ ps\ has 
the form of a standard orientation test in dimension <5). Therefore to complete our analysis 
we must: 

1. Analyze the statistical behavior of \p1p2 ■ ■ .ps\\ 

2. Analyze the statistical behavior of power -(ps+i, S); 

3. Obtain a convenient upper bound to the product of two random variables. 

These tasks are the object of the next three subsections. The main idea of the proof is to use 
the fact that W = power(ps+i, S) does not depend actually on pi,P2 ■ ■ -Ps but only on their 
circumscribing sphere. 



2.1 Orientation test 



In [DP95] we have shown that, given 8 points uniformly distributed in the unit ball Bs in 
dimension S, 

Prob(\pi,p 2 . . .ps\ < V) < osV 

v 5 _ 

where as = 5 {} and Vj denotes the volume of the unit ball in dimension j. 

V 6 

In fact, these results can be extended without any difficulty to the case in which the value 
of \pi,P2 ■ ■ -ps\ is constrained to an interval [V, V + dV], by simply changing in Equation (c) of 
[ DP98f the integration bounds from f mm< - v ' as - 1 ^ to r mln ( v ' + ' iv '' a < 5 -i) This trivial modification 
readily yields 

Prob(V < \ Pl , P 2 . . .p s \ <V + dV I pi,p 2 . . .ps e Bs) < o-sdV (3) 
Th is result generalizes to the uniform distribution in the unit cube Cs ~ [—1,1]* as in 

Prob(V <\ Pl ,p 2 ...ps\<V + dV \ pi,p 2 ...ps eCs) <ipsdV (4) 

6(6-1) 

where ip s = ■ 

2.2 Power of a point with respect to a sphere 

Given a sphere S, with center c and radius r, we wish to compute the probability for a random 
point p to have a small (absolute value) power with respect to S. 
For a small value V we observe that 

(\power(p,S)\ = \\cp\ 2 -r 2 | < V) =>• ( r - < V V 2 - V < \cp\ < Jr 2 + V < r + ¥-) 
v 1 1 ' \ 2r 2r J 



Therefore the value of the power of p with respect to S is smaller than V if p belongs to a 
spherical crown of S of width — . Clearly, the volume of such crown is given by the measure 
(area) of S multiplied by ^, i.e., it is given by Svsr 6 ' 1 ^ — 6vsr s ~ 2 V ( this holds in our 
hypothesis of small V). 

Thus Prob(power(p, S) < V) is bounded as follows: 
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V a a + da 1 a 

Figure 1: Upper bounding event ab <V 
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The term volume(crown n fi) is the product of 2- by the area of S fl O. At this point we 
assume !] C Cs, which is obviously verified when Q is either Bs or Cs- If r < 1 we bound from 
above the volume of the crown by Sv s r s ~ 2 V < 5v s V. If r > 1 we restrict ourselves to the 
portion of the crown internal to Cs and obtain area(5 PI Cs)^r < area(S PI Cs)V < 5v$V. 

In conclusion, we have 

Prob(power(p,S) < V | S given; p £ Bs) < —V = SV (5) 

vs 

Prob(power{p,S) <V\S given; p £ Cs) < ^-V (6) 

2.3 Product of two random variables 

To complete the analysis outlined above, we need a technical result concerning the probability 
of a product of random variables. 

Let a and b be two random variables such that the marginal probability of a satisfies 
Prob(V < a < V + dV) < AdV and the probability of b conditional on a satisfies Prob(b < 
V\a) < BV, for some constants A and B. Notice that our random variables \pi,Pi ■ ■ -Ps\ an d 
power(p, S) fit the specifications of a and b, respectively. We shall bound from above the event 
ab < V by a union of events of the kind a < a < a + da and b < — , as illustrated on Figure |l[ 

Thus we have 



Prob(ab<V) < Prob{a<V)+ I Probla = a)Prob(b < —la = a)da + Probtb < V) 

a 



< (A + B)V+ / ABV 



r da 



< (A + B)V + ABV In — 



(7) 



Notice that for A and B both > 2 and for V < 1/e, the first term is dominated by the 
second one. 



3 Completing the analysis 

In this section, we present the main conclusion of this note. Recalling that 

A s = [pip 2 ■ ■ .p«[.jxH/jer(pff + i,sphere(pij>2 . ..ps)), 
and the previous bounds, we obtain for the two domains: 

1 



Prob(A s <V\pi...ps+i eB s ) < {a s + 5) V + a s SV In 



V 



Prob(A s < V]pi . . .ps+i £ C S ) < ^Vln I + fa + ^ 



r 



(8) 
(9) 



which express a bound nearly linear in V for t he absolute value of the incircle test for 5 > 1. 
For small values of 5 we recall from [DP98] the (approximate) values of v$, as and ipg- 
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These analytical results can be compared with the experimental results mentioned earlier. 
The latter have been obtained using random point selection in Bs, and are shown in Figure 
|i| They confirm the sublinear behavior for 5 = 1 and a basically linear behavior for S > 2 
near V = 0. However, the constants reported above are far from tigh t when the dimension 
increases, which is a clear byproduct of the technique of proof used in [DP98] to bound %f)s. 



4 Example: 

arithmetic 



3D insphere test with double precision floating point 



We now consider a practical implementation of the insphere test in three dimensions. The 
corresponding expression is given below. We assume that entries (point coordinates) are 
floating point numbers in the range [—1,1] and that they are stored as double precision numbers 
with a 53-bit mantissa. We assume that the computation complies with the IEEE 754 norm. 
We first detail the formula for the insphere test: 
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We now estimate the maximum a priori round-off error using the following standard 
rules: error(x + y) < error(a:) + error(y) + (x + y)2 -54 and error(xy) < x x error(j/) + y x 
error(:r) + xy.2~ 54 . Each computation is analyzed in terms of the elementary operations of 
addition/subtraction or multiplication. 
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If the points are uniformly distributed in the unit cube and snap-rounded to the nearest 
representable point, then the above calculations show that if the insphere test gives a result 
larger than 129 2 -51 (in absolute value), then its sign is reliable. 

For simple precision numbers with 24 bits of mantissa, an analogous statement can be 
made for results larger than 129 2 -22 ~ 2 -15 . 

These results enable us to estimate the probability of failure of such filter, i.e., prob(failure) < 
32(V) In i + 22(V) with V = 129. 2 -51 or V = 129. 2 -22 for the two cases. 

Claim: If the absolute value of the insphere test in three dimensions for points in 
the unit cube computed with 53 (resp. 24) bit arithmetic is larger than 129 2~ 51 < 
6 KT 14 (resp. 129 2~ 22 ~ 3 l(T 5 j then the sign is reliable. The probability of 
failure of the certifier is less than 6 10 -11 (resp. 0.011,). 
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